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POINT CHARGES AND POLYGONAL LINKAGES 


GIORGI KHIMSHIASHVILI*, GAIANE PANINA**, DIRK SIERSMAt, VLADIMIR ZOLOTOV* 


Abstract. We investigate the critical points of Coulomb potential of point charges 
placed at the vertices of a planar polygonal linkage. It is shown that, for a collection 
of positive charges on a pentagonal linkage, there is a unique critical point in the set 
of convex configurations which is the point of absolute minimum. This enables us to 
prove that two controlling charges are sufficient to navigate between any two convex 
configurations of a pentagonal linkage. 


1. Introduction 

We deal with the Coulomb potential of a system of point charges placed at the vertices 
of a planar polygonal linkage. The ultimate goal is to establish the possibility of controlling 
the shape of a pentagonal linkage by varying the values of two charges at the vertices. 
Following a paradigm developed in our previous paper [10] we consider Coulomb potential 
as a meromorphic function on the planar moduli space of linkage and investigate its critical 
points. The basic assumption and motivation for accepting such an approach is that a 
vertex-charged linkage subject only to Coulomb interaction of its charged vertices should 
take the shape with the minimal Coulomb potential. 

The approach developed in [10] was suggested by some recent research concerned with 
the control of nanosystems and other systems with several degrees of freedom [8] , [13] . This 
setting suggests several aspects and problems. In the present paper, we concentrate on the 
following scenario. 

Given a convex planar configuration of linkage, we wish to find the vertex charges such 
that the global minimum of the arising Coulomb potential is achieved at the given configura¬ 
tion. Such a collection of charges will be said to stabilize the given configuration. We prove 
that any convex configuration of pentagonal linkage has a stabilizing system of charges which 
depends continuously on the configuration. This allows one to navigate any (initial) convex 
configuration to any other (target) convex configuration along any prescribed path in the 
space of convex configurations. We will refer to this situation by saying that a pentagonal 
linkage admits a complete Coulomb control. 

We wish to add that this research arose as a natural continuation of our previous joint 
results on Morse functions on moduli spaces of polygonal linkages m, Em- 

The paper is organized as follows. 

Section [2] contains necessary preliminaries. Among them is a formula for the partial 
derivatives of Cayley-Menger determinant (Theorem [1]) which seems interesting for its own 
sake and provides one of our main technical tools. Section [3] briefly sketches an outline of 
the proof of the Coulomb control for quadrilateral linkages. In a sense, it is a prologue which 
introduces the approach to be used for pentagonal linkages in the subsequent sections. 


Key words and phrases. Mechanical linkage, configuration space, moduli space, Coulomb potential, 
Cayley-Menger determinant, Coulomb control. 


1 



2 GIORGI KHIMSHIASHVILI*, GAIANE PANINA**, DIRK SIERSMA^, VLADIMIR ZOLOTOV* 

In Section |4] we charge a pentagonal linkage by a five-tuple of positive charges 
q = (gi, ..., 55 ). We prove that for any collection of charges, the Coulomb potential has a 
unique critical point (which is the minimum point) in the space of convex configurations. 

In Section[5]we control a pentagonal linkage by just two positive charges. This means that 
we put charges (gi, 92 , 94 , s), assuming that 91 , 92,94 are some fixed positive charges, and 

that we can vary s and t. We prove that for any convex configuration P there exists a unique 
stabilizing pair of positive charges s, t. This yields our main result (Theorem|4]) stating that 
these two charges provide a complete control on the space of convex configurations. 

To simplify and properly structurize the presentation all technical proofs are placed in 
the separate Section [ 6 ] 
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session in CIRM (Luminy) in January of 2015. The authors acknowledge the hospitality 
and excellent working conditions at CIRM. We also thank Pavel Galashin for useful and 
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2. Notation and preliminaries 

2.1. Polygonal linkages and their moduli spaces. A polygonal linkage L is defined by a 

collection of positive numbers I = (Zi, called sidelengths, which we express by writing 

L = L{1). Physically, a polygonal linkage is a collection of rigid bars of lengths k joined in a 
cycle by revolving joints. It is a flexible mechanism which can admit different shapes, with 
or without self-intersections. 

By M{L) we denote the moduli space of planar configurations, that is, the space of all 
polygons with the prescribed edge lengths factorized by isometries of 

M{L) = {{pi, ...,Pn)\p^ G \piPz+l\ = k, \PnPl\ = ln}/Iso{R^). 

This is not exactly the moduli space M[L) treated in [5] and [5], where the space of 
polygons is factorized by orientation preserving isometries. However there is a two-fold 
covering A4(L) —>■ M(L). 

By M‘^{L) we denote the set of all strictly convex configurations. We exclude here non- 
strictly convex polygons, that is, those having (at least) one angle equal to tt. The latters 
obviously form the boundary i9M‘^(L). 

The set of all convex configurations, that is, the closure of [L) in the ambient space 
M(L), is denoted by = M^{L). 

In this paper we basically deal with n = 4, 5. For a 4-bar polygonal linkage, M{L) is a 
(topological) circle, whereas M^{L) is homeomorphic to a segment. For a 5-bar polygonal 
linkage, M{L) is (generically) a surface, whereas M^{L) is homeomorphic to a disk D^. 

2.2. Partial derivatives of Cayley-Menger determinant. We present now the defi¬ 
nition and certain properties of the Gayley-Menger determinant which is one of our main 
tools. 

Let Ai, A 2 , A 3 , A 4 be four points in K^. Denote the distances between the points and the 
vectors as 


-A- 


dij — [A^Ajj, d ij — AjAj. 


The Gayley-Menger determinant of the quadruple of points is defined by the formula 
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As is well known, D is equal to 288 times the squared volume of the tetrahedron Ai,A 4 . 


Denote by Sijk the oriented area of the triangle AiAjAk, that is, the orientation is taken 
into account. As an example, Sijk = —Sikj. 

We also need the vectors 

-A _ “t 

i-ik — ^ a i-j y. a jkj 


for which we have 


Theorem 1. 


\S ijk\ — \Sijk\- 

(1) For any four points Ai,..., A 4 in . 
1 d 


we have: 


2ddl^ 


D = —16( S 124, S 234 ), 


where .) denotes the standard scalar product in the space 
(2) In the case where the four points are coplanar, the formula reduces to 
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d'li ~ di2 ~ 1^24 


“24 

“ ^^23 “ 

-2di4 

dh - dh - dh 

— 2^24 


= -4 


(^ 12 ,^ 23 ) C^23i~^24) 

C ^ 12 ,~<^ 24 ) ("1^24, "1^24) 


= —4('^12 X ~^24,~^23 X ~^24) = —16(^124,^234), 

where ” x ” denotes the cross-product of vectors and the last line of equalities follows from 
the Binet-Cauchy formula. 

□ 
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This theorem seems to be new. To the best of our knowledge, this formula appeared in 
the literature only for the case when the points Ai,Ai are coplanar. For example, in [4] 
it is given without proof but with the following comments. ’’This is an important formula 
concerning the Cayley - Menger determinant. This formula is only valid when restricting 
to planar configurations. The minus sign in equation is not included in Dziobek’s original 
paper [5] nor in several later works that utilize the Cayley-Menger determinant. However, 
checking equation on the square configuration indicates the need for the minus sign. The 
correct formula appears in the doctoral thesis of Hampton [7]”. 

We actually use it only in the coplanar case for describing relations between diagonals 
of polygons. However we presented the proof of a more general statement for completeness 
and convenience of the reader more so that [7] does not seem to be easily accessible. 

2.3. Coulomb potential. Placing a collection of point charges qi at the vertices A^ of 
a configuration and considering the Coulomb potential of these charges we get a function 
defined on M{L). We will refer to this setting by speaking of a vertex-charged linkage with 
the system of charges q= (gi,..., q„). 

Recall that the Coulomb potential E oi a. system of point charges € K placed at the 
points Pi of Euclidean plane is defined as 

E = = ( 1 ) 

... 0,ij 

where dij = is the distance between the i-th and j-th points. 

Since we are only interested in critical points of Coulomb potential, addition of a constant 
makes no difference. By the very definition of polygonal linkage, the distances corresponding 
to two consequent vertices in formula (ED remain the same for all conhgurations of linkage. 
Hence their sum is constant for any fixed collection of charges for our purposes it is sufficient 
to work with the effective Coulomb potential E defined as 

E = ( 2 ) 

i<j-i 

where dij is the length of diagonal between (non-neighboring) i-th and j-th vertices of the 
configuration. 

We say that a collection of charges stabilizes the configuration P A E attains at P its 
minimal value. In this case we say that P is the minimum point of E. 

We explicate now the setting and notation in the cases considered in the sequel. For 
n = 4, we put one positive charge t at the first vertex. The rest three vertex-charges are -1-1. 

For n = 5, we put two positive charges s and t at any two non-neighboring vertices and 
say that s and t are controlling charges. The rest three charges are some fixed positive 
numbers. 

It is proven in [10] that, for positive charges s, t, the global minimum of E on the moduli 
space M{L) belongs to . 

3. Coulomb control of quadrilaterals 

As a visual illustration of the paradigm we briefly recall the main results of [TU]. For a 
quadrilateral linkage, we put charges equal to 1 at three vertices, whereas the fourth vertex 
is charged by t. Then we have 



t 

^24 
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The following facts were established in [TO] : 

(1) For a given convex quadrilateral P £ M{L), there exists a unique t such that P is 
a critical point for this charge. In this case, t is positive. 

(2) For a quadrilateral linkage, E has a unique minimum in . 

This means that, for a quadrilateral linkage, we have a complete control of convex con¬ 
figurations by point charges at its vertices. 

The ’’picture of what happens” which is behind the proofs of [10] is such. The configu¬ 
ration space of a quadrilateral is a circle, and therefore the relation between the diagonals 
is important since both diagonals appear in E. We remind that four points Ai ,..., A 4 are 
coplanar whenever their Cayley-Menger determinant vanishes. Therefore the squared diag¬ 
onals X = d\^,y = ^24 of the quadrilateral are related by a cubic equation D{x, y) = 0. The 
equation defines a real elliptic curve lying in the {x, y) plane. It contains two connected 
components: one is unbounded and the other one is a convex closed curve, which we denote 
by C. We make use of the mapping 

'H -.M 

«'(P) = {x,y) = (di3,d^4), 

which maps a configuration of the linkage to the point in the plane whose coordinates are 
squared diagonals. The configuration space is mapped to the bounded component of the 
elliptic curve, whereas the space of convex configuration is bijectively mapped to some arc 
of the curve C. 

The rest of the proof from [TOj uses the idea of convexity. Namely, we can consider E as 
a convex function defined on the (x, y)-plane. We are interested in its restriction to the arc 
of the convex curve C. The restriction of a convex function to a convex curve is (in general 
case) not necessarily convex. But in our case, due to the specific position of the gradient 
VP with respect to the mentioned arc, and with a special parametrization, the restriction 
is a convex function, and therefore has a unique critical point which is the minimum. 

In Section |4] we apply a similar approach: we introduce a mapping $ from the configura¬ 
tion space of a pentagonal linkage to the space by mapping a configuration to its squared 
diagonals. The image is a 2-dimensional manifold with corners, which appears to be the 
intersection of convex hypersurfaces. 


4. Critical points of Coulomb potential 


Now we assume that we are given a polygonal 5-bar linkage L with the vertices denoted 
by Ai,...,A^, as in the previous section. We also adopt notations for and bi from 
the previous section. Since it is a linkage, we now assume that the sidelengths IA 1 A 2 I = 
tti,..., 1 ^ 5 All = 05 are fixed. Let 


E 

i<j-i 


\AiAj 


be the (effective) Coulomb potential of the positively charged linkage L. Here we denote 

Cij .— QiQj ^ 0 . 

Let M^{L) be the set of all convex configurations of L. The main result of this section 
establishes the uniqueness of critical point of E in M^{L). 


Theorem 2. The potential E has a unique critical point in the set which is the 

(global) minimum of E on M{L). □ 


The rest of the section is devoted to the proof of this theorem. 
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Figure 1. The space M^{L), the polar curve T, and its components Ti and r2. 


Definition 1. A k-slice C M^(L) with respect to the diagonal ^3^5 is the set of all 
(convex) configurations such that 64 = IA 3 A 5 I = k, where k is some constant. 

Each slice is an analytic curve homeomorphic to either a line segment or a point. There 
are exactly two slices, called terminals which are points. The (disjoint) union of all the slices 
equals the set M^{L). 

Proposition 1. For each of the slices Xk, there exist two possibilities: 

(1) The restriction of the potential E\xk has a unique critical point which is the mini¬ 
mum point. 

(2) The restriction of the potential E\x^. has no critical points. Then the minimum 
value is achieved at the boundary of the slice. 

In any case, the minimum point of the restricted potential is unique and depends continuously 
on k. □ 


Now we pass from one particular slice to the set of all convex configurations. 

Definition 2. On each of the slices, we mark the minimum point of E\xk- Taken together 
for all slices the marked points form the polar curve T C M^{L) of the potential E with 
respect to the diagonal A^A^. 

The polar curve is a piecewise analytic curve homeomorphic to a line segment. Its end¬ 
points are the terminals. The polar curve may also contain some segments lying on the 
boundary of M^{L). The intersection of the polar curve with the interior part M^{L) is a 
finite set of connected components T^. 

Since all critical points of E belong to the polar curve, to prove Theorem it suffices to 
show that E\r has a single critical point in T. 

Proposition 2. The potential i?|r has a unique critical point in each of the connected 
components T^ of the polar curve T. □ 


With these ingredients at hand, we are ready to prove the main result. 
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Proof of Theorem [ 2 l Let us start with the most symmetric case: the equilateral linkage 
and equal charges. We prove that E has a unique critical point in M^{L). 

Indeed, each slice Xk contains an axial symmetric configuration, that is, there exists a 
convex polygon 711 ^ 2 ^ 3714^5 such that: 

= k, 


ZA2A3A4 — /1A4 A^Ai, AA1A2A2 — AA^AiA2 . 

By symmetry reasons, A1...A5 is a critical point, and therefore, (the unique) minimum 
point of the potential E restricted to Xk- If Xk is not a terminal point, the critical point 
lies in the relative interior of Xk- Therefore, for this particular case, the polar curve lies in 
the space and by Proposition [2] we obtain the result for the equilateral pentagon 

and equal charges. 

For an arbitrary pentagonal linkage L and arbitrary (positive) charges, we use reductio ad 
absurdum. Assume that F has more than one component, and E has more than one critical 
points. We start continuously deforming the lengths and the charges aiming to the above 
symmetric linkage. Let us inspect the behavior of the critical points during the deformation. 
During this process all objects change continuously: the spaces M{L) and the slices 

Xk, the curve F and the critical points. The number of components of F can increase, but 
also decrease (as soon as the polar curve crosses the boundary). In the beginning of the 
deformation we have more than one critical points, whereas at the end the critical point 
is unique. The confluence of critical points cannot happen in the interior, since then two 
branches of F have to be connected, but the limits of the critical points are distinct on that 
new branch unless they meet at the boundary. However two critical points on one branch is 
impossible by Proposition [2j 

Moreover, it is proven in [10] that E has no critical points on the boundary of M^{L). 
This means that change of the number of critical points is not possible, so from the beginning 
we have only one critical point which is the minimum. □ 

5. Coulomb control of pentagonal linkages 

For a generic pentagonal linkage, we put charges s and t to the vertices 5 and 3 respec¬ 
tively. For this case we have 

E = I I g25 ^ st 

dlA ^24 dl3 ^25 ^35 

Now we wish to understand whether two charges can provide a complete control of convex 
pentagons. We begin with presenting a simple but conceptually important general obser¬ 
vation valid for arbitrary polygonal linkages. This observation makes essential use of the 
specific form of Coulomb interaction and underlies much of the further discussion. 

Proposition 3 . For any n-gonal linkage L and any configuration P € M(L), the stabilizing 
charges for P are solutions to a system of n — 3 quadratic equations in n unknowns with the 
coefficients algebraically expressible through the lengths dij of the diagonals of P. 

Proof. The proof is obtained by a standard use of the Lagrange multipliers method 
for constrained optimization of £1 as a function of diagonals dij. In our case the number 
of variables dij is equal to k{n) = n{n — 3)/2. It is also easy to see that the number of 
independent constraints Dj = 0 given by the Cayley-Menger relations for the diagonals is 
equal to l{n) = (n —2)(n —3)/2. Following prescriptions of the Lagrange method we consider 
a {l{n) -I- 1) X fc(n)) functional matrix J{P) having the gradient of target function Eq as the 
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first row, and gradients of constraints Dj = 0 as the remaining rows. Notice that the values 
of charges qj appear only in the first row. 

According to Lagrange criterion a system of charges q yields a constrained critical point 
of E on M(L) if and only if the rank of J{P) is not maximal. In other words, all 
{l{n) + 1) X {l{n) + l)-minors of J{P) should vanish at point P which gives us a system 
of algebraic equations S{n) for q each of which is of degree not exceeding two. By linear 
algebra the number of independent equations in the system S{n) is equal to n — 3. The 
statement about coefhcients can be verified directly. 

□ 

For a given P S M{L), let us denote by Stab{P) the set of solutions to S{n), i.e. the set 
of all charges stabilizing P. For obvious reasons, one may await that generically Stab{P) is 
three-dimensional and it is easy to see that this is true for n = 4. Indeed, in this case we have 
just one homogeneous quadratic equation in four unknowns qi of the form Agiga -|- Bq 2 q 4 . 
It follows that the solution set Stab{P) is a cone over the one-sheeted two-dimensional 
hyperboloid. We add that similar results can be obtained for n > 5 using a general method 
for geometric and topological investigation of intersections of real quadrics developed in [ 1 ] . 

In general, the structure of Stab{P) can vary from point to point, which makes it unclear 
how to navigate from one configuration to another one. Developing effective navigation 
would be easier if the sets of available charges were finite. For dimensional reasons the 
number of solutions to S(ri) may be finite if we reduce the number of unknowns (charges) 
to n — 3. So it becomes natural to suggest that Coulomb control should be possible if the 
number of controlling charges is n — 3, which will be assumed from now. In this setting 
Proposition[3]yields a modified system S'{n) and a universal upper estimate for the number 
of stabilizing charges which follows merely from Bezout theorem. 

Lemma 1. If Stab{P) is finite then the number of stabilizing charges for P does not exceed 

2n-3_ □ 

Examining the structure of system S'{5) we conclude that if Stab{P) is finite then its 
cardinality in fact does not exceed two. The reason is that if the controlling charges are not 
neighboring then the system S"(5) consists of one linear and one quadratic equation. In this 
case explicit expressions for the coefficients can be found in Section (|6.4I) . 

Moreover, if we only consider convex pentagons with two controlling charges then it turns 
out that there exists exactly one pair of positive charges stabilizing P, which is the desired 
situation for our purposes. To prove the latter fact we need one more lemma. 

Lemma 2. [10] For st < 0, a convex pentagon is never a critical point of E. □ 

Theorem 3. For each convex configuration P of a pentagonal linkage, there exists exactly 
one {s,t) G (M>o)^ such that P is a critical point for the charges {s,t). □ 

We are now able to explicate the complete control for convex configurations of arbitrary 
pentagonal linkage, which is the main conceptual result of this paper. 

We describe first our system in the general context of state space and control space. Our 
state space is M{L), the space of pentagons with length vector L. The control space is the 
subspace of charges {qi, q 2 ,t, 94 , s). So we work with as control space. Given a potential 
function 

E = Eg : M(L) X ^ R, 

which maps {P,q) to Eq{P), one can consider the critical surface C := {{P,q)\VpE = 0}. 
We also consider the subset of minima Cmin- 

The restriction to C of the projection to each of the two factors of M{L) x gives us 
two maps: tti : C —)• M{L) and 712 : C —>• R^. Generically 712 is a finite cover which is a 
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Figure 2. The critical space C, the state space and the control space. 

local diffeomorphism away from the bifurcation set in K^. A main question is now if it is 
possible to connect the two configuration Pq and Pi by the lift (via tti) of a path in the 
control space; if possible through the space Cmin- If so, it is important to construct such a 
path, which we could call an explicit control to navigate from Pq to Pi. 

Theorem [2] about the uniqueness of the absolute minimum gives an affirmative answer if 
Pq and Pi are convex. Indeed ( 712 )“^ o tti is a bijection between {L) and K?>o- 

The corresponding algorithm to connect Pq and Pi can be described more precisely as 
follows. First, we calculate the starting and target stabilizing charges (so,fo) and (si,ti) 
using formulae coming from the proof of Theorem 16.41 in Section 6.4. Next, we choose a 
path joining (so,to) and (si,ti) in the control space R^>o- So if we continuously change the 
values of stabilizing charges along the chosen path the configuration will move to the target 
configuration through convex configurations without meeting a bifurcation point. 

We summarize the above as follows. 

Theorem 4. Assume that a pentagonal linkage L is charged by (gi, 92 , t, 94 , s), where 
< 71 , 92,94 are some fixed positive charges. For any starting convex configuration Pq G M^{L) 
together with any target convex configuration Pi G M^(L), there exists an explicit control 
by two positive ruling charges {s,t) which yields a path between Pq and Pi lying in M^{L). 

Remark 1. Since there exist infinitely many paths joining the starting and target stabiliz¬ 
ing charges in the parameter space, several natural problems may now be formulated and 
explored in our setting. In particular, one can consider various problems in the spirit of 
optimal control. For example, since there exist various natural Riemannian metrics on the 
configuration space one can investigate which connecting path in parameter space gives the 
shortest path in the configuration space. Thinking of linkage as a device performing certain 
task one may wish to specify a path in the working space of a certain vertex, say, to avoid 
collision with an obstacle. Further examples of such problems can be easily formulated and 
will be discussed elsewhere. 

Remark 2. It seems worthy of noting that placing the controlling charges at adjacent ver¬ 
tices the situation becomes worse. Namely, we will not be able to reach all of the convex 
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polygons (see |10j ). For instance, we will never be able to have two vertices simultaneously 
aligned for an equilateral pentagon. By continuity reasons an entire neighborhood of such a 
polygon becomes unreachable. 


6 . Proofs 

6.1. Notation and applications of Cayley-Menger determinant’s derivatives. We 

start with an elementary application of Theorem [T] about Cayley-Menger determinant. 

Assume that we are given five points Ai,...,A^ in the space Let us introduce some 
ad hoc notation convenient for our purposes: 

bi ■= \Ai_iAij^i\, 

Xi = b1, 

where we treat the indices modulo five. 

We think of at as the edges and of bi as diagonals of the pentagon Ai,...,A^. Notice that 
we permit that the pentagon may be non-planar and non-convex. 

We also need the squared diagonals Xi which will play the role of variables. 

Let Di be the Cayley-Menger determinant for the vertices of quadrilateral which is ob¬ 
tained by cutting of a triangle along i-th diagonal. For example, is the Cayley-Menger 
determinant for the quadruple Ai, A 2 , A 3 , A^: 



0 

1 

1 

1 

1 


1 

0 

aj 

X2 

4 

Di = 

1 

al 

0 

4 

Xi 


1 

X2 

4 

0 

X 4 


1 

4 

Xi 

X 4 
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Applying Theorem [T] in this situation we get several useful equalities. 
Lemma 3. Assume that Ai ,..., A 5 are coplanar points. Then we have: 


- Da = — 16S'i25<S'235- 

2 0 x 2 

1 

— Da = +16S'i25<S'i23- 
2 0 x 4 

77-5- D4 = —16S'i23S'i35. □ 

2 axi 


6.2. Proof of Proposition [H Throughout the paragraph we fix 64 and parameterize the 
slice Xk by X 2 . In this setting the lengths of diagonals bi, the squared lengths of diagonals Xi, 
and the restriction of the potential E\xi, are the functions in the variable X 2 . All derivatives 
(denoted by’’prime”) mean derivatives with respect to X 2 . For instance, we write (&i)' = 


Lemma 4. (1) 


{-xi)' 


<S'i25 <S'235 

5' i 23 <5135 


^234 5i25 

5i35 5i23 


(2) 


(-a’ 3 )' 
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(3) 


{xbY = 


'S'145 


’135 


Proof follows from Lemma |3] by using the implicit differentiation formula. 


□ 


Lemma 5. For a fixed slice Xk, we have (bi)' < 0, ( 63 )' < 0, ( 65 )' > 0 on the relative 
interior of Xf^. This implies that bi J,, 63 f, 65 f- 

Proof follows from Lemma IH □ 

Lemma 6 . For a fixed slice Xk, we have 

(1) -(xi)" > 0 

(2) -{x,)" > 0 

(3) gfft 

Proof. (1) follows from second order relation between the diagonals of quadrilateral, see [lO] . 
The statement (2) can be proven the same way as (1). 

By Lemma m 

{xsY _ >§'234 
(a;i)' -S '235 

The statement (3) follows now from direct computation of the derivative. □ 

Lemma 7. For a fixed slice Xk, we have >0 on the relative interior of Xk- 

Proof. The Coulomb potential for the pentagonal linkage writes as: 

E = C 25&1 ^ + C 13&2 ^ + C 24&3 ^ + C 3564 ^ + C 1465 
The derivative of the potential along the slice Xk writes as: 

E\'xk = C25bf^i-XiY + C2ib^^{-X3Y - Cisbf^ - Ciib^^ixsY- 
Using the fact that (xi)' < 0, we can rewrite this formula as: 

^\'xk = {-XlY {c25K^ + C2 ib3^ (-Xi)' ) ~ - C14&^^(X5)'. 

To prove that E" > 0 it suffices to show that: 

(-xi)' > 0 , 

(c25K^ + C24^3 ^ ^ 

(-^ 1 )" > 0, 

(c256r^+ 02463 -" 1 ^^) t, 

(- 01362 ^) t) 

(-Ci 46^^(X5)') t ■ 

All these statements we know from Lemma [5] and Lemma [ 6 ] □ 


Lemma[7]implies Proposition[T]straightforwardly since E is strictly convex (as the function 
in the variable X 2 ) on each of the slices. 
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6.3. Proof of Proposition [2l As was already mentioned, M^{L) is a two-dimensional 
closed disk. We embed M^{L) in by mapping each configuration to the squared lengths 
of its diagonals Xi,... 

$ : M{L) 

This mapping sends M'"{L) to its image bijectively so this is indeed an embedding. The 
mapping extends by the same rule to the entire configuration space. However on the entire 
configuration space it is not injective. 

Now we think of i? as of a function defined on R^. 

K — £li 

~ bi b 2 bs bi bs' 

We will deal with the signs of components of its gradient. As a matter of fact they all 
are negative: 


VF — ^13 C 24 C35 C41 N 

2\(bl)i' {bl)i' {bl)i\bl)r ^ ; ; ; ; J. 

Here we denoted by the set R_ x R_ x R_ x R_ x R_. In the sequel we 

use analogous notations regarding various combinations of signs of expressions in question. 

Let 7(t) = (7i(t),72(i),73(i),74(i),75(i)), where t G [0,1] be a C^-smooth curve in R®. 
Later we shall assume that 7 is the polar curve but now we consider just any smooth curve. 
We have 


Ehit)) = 




( 7 fc)- 


where the sum is over all triples such that i -\- 1 = k(rnod5), fc -|- 1 
by prime ' the derivative ^ and compute 


j{mod5). We denote 


E{i{t)y = -\Y. 


Cjji'Jk) 

ilk)^ 


E{l{t))" 



c. 4 -[( 7 fc)f 

(7fe)5 


From this we conclude that if the function 


1 ’^iji'lk)" 

2 ilk)^ 


1 

2 (7fe)5 


(vl;, ( 7 )") 


is non-negative on the curve 7, then E has a single critical point on the curve 7. 

We remind that we denote by Di the Cayley-Menger determinant for the points A2, A3, A4, A5, 
denote by the Cayley-Menger determinant for for the points Ai, A3, A4, A5, and so on. 

The system 

Di = D 2 = ... = = 0 

defines a surface which contains the image of M'^[L) under the above described mapping. 

Let us consider the function D2 separately. It is a polynomial of degree 3 in variables 
X2TXi,x^. It doesn’t depend on a:i,a:3. 


Consider an open 3-arm R with edgelengths 03,04,05. It is a subchain of our 5-linkage. 
The configuration space M{R) of the arm is homeomorphic to the 2-torus. 

We map the configuration space of the arm to R^ by mapping each configuration to the 
squared lengths of its three diagonals 64,65 and 62: 


$2 : M{R) R3. 
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The image of $2 belongs to the set D 2 = 0. 

We denote by [R) those configurations of the arm that are subconfigurations of some 
element of , that is, that are extendable to a convex pentagon. 

Lemma 8. (1) $2 maps [M'^{R)) to its image bijectively. 

(2) $ 2 (Af^(i?)) is a convex surface (with boundary) in 

(3) The map $2 on M^{R) has maximal rank except for aligned configurations. 

Proof. The surface ^ 2 {M{R)) is a closed surface contained in D 2 — 0. The surface 
$ 2 (M(i?)) bounds in some body. The image of $2 belongs to the set D 2 = 0. D 2 is a 
polynomial of degree three, therefore each generic line intersects 132 = 0 at at most three 
points. Since a line intersects a closed surface at an even number of points, each generic line 
intersects ^ 2 {M{R)) at most at two points. □ 

Notice that this surface is contained in a ’’box” in the positive octant. 

From Theorem [T] we know the signs of all entries of all the D^. For example, 

VD 2 = (0; 16S'345«S'i 45; 0; —16S'i34S'i45; —16S'345«S'i35). 

We present all these signs in the following table: 


VZ3i 

+ 

0 

- 

- 

0 

VZ32 

0 

+ 

0 

- 

- 

VZ33 

- 

0 

+ 

0 

- 

VZ34 

- 

- 

0 

+ 

0 

VZ35 

0 

- 

- 

0 

+ 


Lemma 9. The gradient VI32 is the inner normal vector of the surface D 2 = 0. Similar 
statements hold for other Di. 

Proof. It is sufficient to check if the gradient points inside or outside for one point only. 
We assume that we pick a pentagon without any aligned edges. We reduce the dimension 
as follows. First fix X 2 and after that X 4 . There are only two quadrilaterals satisfying this 
condition: one non-convex and one convex (which has bigger X 5 ). The intersection with 
the convex body is an interval. The x^ component of the gradient is there negative, so the 
gradient vector points inside. □ 

Remark 3. The curve F n [L) is given by E' = 0 and therefore is an algebraic curve. 
(We remind that E' is the derivative of E along the slice parameterized by X 2 ). Lemma 
claims that E" > Q on X^. Hence, we conclude (with the implicit function theorem) that 
T r\ {L) is a smooth curve which intersects X}~ transversally. 

Let 7 ( 3 ) = 7 (x 4 ) parameterize the image $(F) of the polar curve. 

We are now able to prove Proposition [2] by establishing the inequality 

( VE , 7") > 0 on each of the branches Fj. 

Assume that a point 7(1) does not lie on the boundary of The curve 7 lies on 

$(M‘^(L)) which is the part of the convex component of 13^ = 0 for each i. As we explained 
above, this component of 13^ = 0 is the boundary of some convex body, so by the Darboux 
formulae for the normal curvature we have: 


( 7 (t)";VA) >0 
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Let now ^ be a tangent vector of the slice at the point 7(t). By Lemma [ 5 ] we can assume 
that 

C s [~; +; O; +]■ 

By the definition of polar curve we have (Vif, C) = 0, and obviously {\ 7 Di,() = 0- To 
show that (Vif;7(t)") > 0, it suffices to show that (Vif;7(t)" + \() > 0 for some A. For 
any A, we have already shown that {\ 7 Di;j{t)" + X() > 0. 

Since 'y{t) is parameterized by x^, we have 7(t)" G [*; *; *;0; *], where * denotes entries 
of unknown signs. Since (ViA4;7(t)") > 0, only three cases are possible: 

7(t)" G 0;*], 7(t)" G 0; *], and 7(t)" G 0; *] 

We treat these cases separately. 

(1) The first case is simple: since > 0, (VI?2;7(t)") > 0 we have 7(t)" G 

;0;—), and lyE\^{t)") > 0, since VE G This completes 

the proof of Proposition [ 2 ] 

(2) In the second case we use {VD 2 ',j{t)") > 0 and get 7 (t) G (+; —; =i=;0; —). Here we 

have two cases: 7 (t)" G (+; —; +; 0; —) and j{t)" G (+; —; —; 0; —). Assume we have 
7 (t)" G (+; —; +; 0; —) (the other case is treated similarly). Let us take 7 (t)" + 
and look how the signs change when we continuously increase A from 0 to +oo. 
We start from (+;—;+;0;—) and go to C G ;0;+), which means that all 

the entries (except for 0) change their signs. Let us enumerate the signs this way: 
(+i; — 2 ; +3; 0; —5). The inequality {'VD 4 ] 'y{t)" + X() > 0 implies that the first sign 
changes before the second. The inequality (ViAi; 7 (t)" + AC) > 0 implies that the 
third sign changes before the first. {VI? 2 ; 7 (t)” + AC) > 0 implies that the second 
sign changes before the fifth. So at some moment we necessarily have 

7(t)" + AC G (—; —; —; 0; —). Then (VA; 7(t)" + AC) > 0, which implies 
(V£’;7(<)") > 0, and we are done. 

( 3 ) The third case is treated similarly to the second one. 

The proof of Proposition [ 2 ] is now completed. □ 

6 . 4 . Proof of Theorem [ 3 l The equilateral case was already proven in m- We follow its 
proof. We rewrite potential in the form 

E = g2S st 

h b 3 ^ h h bi 

Take now the diagonals 64 and 62 as local coordinates in a neighborhood of P. 

The polygon P is a critical point of E means that dE vanishes: 

dE qiQiai g 294 / 3 i sg 27 i st 

dbi ~ bl + 52 ^ bl ^ hl~ ^ 

and 

dE _ 91540:2 9294/32 qit 92S72 _ „ 

db 2 ^ hi + bl ^ bl^ bl ~ ' 

where 

oi = 965/964, / 3 i = 963/964, 7i = 961/964, and 


02 = 965/962, (32 = 963/962, 72 = 961/962 . 
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We get a system in two variables s and t which reduces to the following quadratic equation 
in s: 

A + Bs + Cs^ = 0 

with 

qiQiai , < 7294/31 


A = 






„ 9271 &2/94Q<2 , 9294/32 s 

^ - 7/2 1^3- + „ la 


b\ hi hi 


qibl 


and C = — 


^27292 


qiblbl ' 

Since AC is negative by Lemma lU the equation has exactly one real positive solution s. 
By Lemma O t is also positive, which completes the proof. 


7. Concluding remarks 

By Proposition |3] for any n-gon we get a system S{n) of quadratic equations for the 
stabilizing charges. The existence and structure of real solutions to this system can be 
analyzed using topological methods of real algebraic geometry [T]. Finding the number of 
positive solutions to the reduced system S'{n) is also possible using methods of [5]. So it 
is still unclear whether a similar Coulomb control is possible for bigger number of edges. 
However, our expectations are: (1) the Coulomb potential has a unique critical point (which 
is the global minimum) in the domain of convex configurations, (2) for a complete control, 
the non-ruling charges should not be put at three consecutive vertices, and (3) a convex 
configuration may have several collections of positive stabilizing charges. 

Theorem [3] implies that we can navigate from any convex configuration to another convex 
configuration along any path joining their stabilizing charges in the space of charges. Since 
the space of charges is convex, one can use just the segment joining the stabilizing charges. 
It will be interesting to visualize the arising movement of linkage in the configuration space. 
This enables one to consider several natural versions of the optimal control problem for 
vertex-charged pentagonal linkages in various contexts. 
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